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D ■ Abstract 

We construct stationary non-equilibrium black funnels locally asymptotic to global 
AdS4 in vacuum Einstein-Hilbert gravity with negative cosmological constant. These 
, are non-compactly-generated black holes in which a single connected bulk horizon 

\ extends to meet the conformal boundary. Thus the induced (conformal) boundary 

' metric has smooth horizons as well. In our examples, the boundary spacetime contains 

. a pair of black holes connected through the bulk by a tubular bulk horizon. Taking 

\ one boundary black hole to be hotter than the other (AT ^ 0) prohibits equilibrium. 

The result is a so-called flowing funnel, a stationary bulk black hole with a non-Killing 
horizon that may be said to transport heat toward the cooler boundary black hole. 
While generators of the bulk future horizon evolve toward zero expansion in the far 
future, they begin at finite affine parameter with infinite expansion on a singular past 
rS ' horizon characterized by power-law divergences with universal exponents. We explore 



(N 



. both the horizon generators and the boundary stress tensor in detail. While most of 

our results are numerical, a semi-analytic fiuid/gravity description can be obtained by 
passing to a one-parameter generalization of the above boundary conditions. The new 
parameter detunes the temperatures Tbuik bh and Tbndy bh of the bulk and boundary 
black holes, and we may then take a = ^"'^^ and AT small to control the accuracy of 

-'bulk BH 

the fluid-gravity approximation. In the small a, AT regime we find excellent agreement 
with our numerical solutions. For our cases the agreement also remains quite good even 
for a ~ 0.8. In terms of a dual CFT, our a = 1 solutions describe heat transport via 
a large N version of Hawking radiation through a deconfined plasma that couples 
efficiently to both boundary black holes. 
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1 Introduction 

We focus here on the classic problem of heat transport far from equilibrium, and away from 
the perturbative regime. If the system of interest is an appropriate strongly coupled large 
N conformal field theory (CFT), we may use gauge/gravity duality to exploit a perhaps- 
more-tractable description as a semi-classical bulk gravitational system. We will consider 
the classical limit in cases where the bulk description may be truncated to A < Einstein- 
Hilbert gravity. Our work complements perturbative computations of heat transport in this 
regime (e.g. [T]), as well as non-perturbative studies of thermalization (see e.g. [21 El II] for 
recent examples and further references) and holographic Shockwaves [5l [6] . 

Suppose in particular that we couple a CFT in d spacetime dimensions to heat sources 
or sinks of finite size and at finite locations. A convenient way to introduce such sources 
is to place the CFT on a background non-dynamical spacetime containing stationary black 
holes with surface gravity k, which have temperatures k/27t due to the Hawking effect. As 
we review in section [2] below, this problem may also be generalized so that the field theory 
temperature at the black hole horizon differs from k/27i. But since no information can flow 
outward across the horizon, the choice of a black hole metric is nevertheless useful to decouple 
our CFT from the details of the heat sources and sinks. The problem of heat transport then 
becomes one of computing the expectation value of the stress tensor in the given background 
with the stated boundary conditions. Since the background spacetime is not dynamical, we 
can choose the metric at will. In particular, we can include as many black holes as we like 
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at locations of our choosing, and we are free to assign their surface gravities as desired. Of 
course, since we consider CFTs, we may also conformally rescale the background metric to 
reinterpret our heat sources/sinks as being infinitely large and located at infinite distance; 
more will be said about this alternate interpretation in section |2] below. 

Gauge/gravity duality for large N field theories has been used to study related settings 
in[8l[a[IOl[IIl[T2l[l3l[Ta[l5l[T6l[I71[l8l[Tg 

In this context, the d- dimensional black hole spacetime on which the CFT lives becomes 
the conformal boundary of a {D = d + l)-dimensional asymptotically locally anti-de Sitter 
(AlAdS) spacetime and we henceforth refer to our heat sources and sinks as boundary black 
holes. Though the above explorations in gauge/gravity duality involved certain tensions and 
subtleties, the picture that emerged in |22j (building on [19j) is one with two important 
phases for each boundary black hole, even when the CFT state is assumed to contain a 
deconfined plasma. See [30] for a condensed review. In the so-called "funnel phase" a given 
boundary black hole is connected to distant regions of the boundary by a bulk horizon along 
which heat may be said to flow (say, if unequal temperatures are fixed at the two ends). 
But there is no such connection in the contrasting "droplet phase." Figure [1] depicts both 
phases for a simple case in which the boundary spacetime is asymptotically fiat. In the 
CFT description, the funnel phase allows a given boundary black hole to exchange heat with 
distant regions much as in a free theory with a similar number of degrees of freedom. One 
may say that grey body factors are 0{1) even at large A^. But in the droplet phase there is 
no conduction of heat between a given black hole and the region far away at leading order in 
large A^. In effect, all grey body factors associated with the black hole vanish at this ordei0, 
so that the black hole does not couple efficiently to the surrounding plasma. Additional 
phases are also possible that conduct heat between subsets of nearby black holes but not to 
infinity. 

Until recently, both funnel and droplet solutions were largely conjectural. Explicit ex- 
amples were known only in rather contrived settings or in low enough dimensions that all 
properties were determined by conformal invariance. But numerical methods were used to 
construct more natural droplets in and more natural funnels in [21] . An interesting detail 
is that the droplet solutions of [28] contain deformed planar black holes (see figure [T]) with 
vanishing temperature. Constructing black droplet solutions that include finite-temperature 
(deformed) planar black holes remains an open technical challenge, though perturbative 
arguments give strong indications that they exist. 

The above funnel and droplet papers largely focussed on cases without heat fiow; i.e., ei- 
ther on droplets (in which heat does not fiow in the approximation that the bulk is classical) 
or on equilibrium funnels. The one exception was [30] which showed that, by changing con- 
formal frames, rotating BTZ black holes [32| 133] in AdSs can be re-interpreted as describing 
heat transport in 1+1 CFTs. Here the standard left- and right-moving temperatures Ti,Tji 
of the BTZ black hole correspond directly to the temperatures of the left- and right-moving 
components of the CFT. Due to the strong constraints of conformal symmetry in low di- 
mensions these components do not interact and the temperatures Tl,Tr must be constants 

""^To be more precise, the grey body factors are non-zero only for a number of degrees of freedom that 
scales like N'^ = 1. 
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Figure 1: A sketch of the relevant solutions: (a) a black funnel and (b) a black droplet 
above a deformed planar black hole. For simplicity, we take both solutions to asymptote in 
the horizontal direction to the so-called planar AdS-Schwarzschild black hole. As a result, 
both describe possible states of a CFT on an asymptotically flat black hole spacetime filled 
with a deconfined plasma at constant temperature. In each figure, the top line corresponds 
to the spacetime on which the CFT lives; i.e., to the conformal boundary of the AlAdS bulk. 
The dots denote horizons of the boundary black holes. The shading marks regions inside the 
bulk horizons. 



if the heat flux is stationary. In addition, the flow of heat is necessarily isentropic (having 
no local generation of entropy). 

We refer to black funnels transporting heat as "flowing funnels." Since none of the above 
special properties should hold for d > 2 {D > 3), higher dimensional flowing funnels should 
be quite different than those found in [30j. For example, a bulk horizon connecting two 
boundary black holes of different temperatures should (at least in some rough sense) be 
describable as having a temperature that varies along the horizon. But recall that there 
is no gwierally accepted definition of horizon temperature which allows this temperature 
to varjo. Indeed, the fact that any definition of temperature should vary implies that the 
horizon is not Killing, which is already a novel property for a stationary black holfl This 
suggests that the horizon generators have positive expansion (though of course tending to 
zero in the far future), so that they caustic at finite affine parameter in the past. It is natural 
to expect this caustic to occur at a singular past horizon [22], and section |5] confirms this 
picture for our solutions. 

We focus below on what we call D = A global flowing funnels, by which we mean de- 
formations of the global AdS4 black string (also known as the Banados-Teitelboim-Zanelli 
(BTZ) black string; see e.g. |10] where the solution was obtained as a special case of the AdS 

^Except of course within the domain of the gradient expansion, as in the fluid-gravity correspondence 
[21]; see also [3S]. For proposals in more general contexts see e.g. [35] for a recent paper and references. Our 
solutions may therefore provide interesting testbeds for such proposals. 

^For compactly generated horizons, this behavior is forbidden by the rigidity theorems [371 |3H1 ISH] ■ But 
our bulk horizon is non-compact since it extends to the conformal boundary. 
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Figure 2: A sketch of a t = const, slice of the BTZ string (11. ip . The two pieces of the 
boundary at z/i^ = ±7r/2 are conformal to two copies of the BTZ black hole, sketched 
above as the two hemispheres of an S^. These boundary black holes are joined at infinity 
(the dashed line around the equator of the sphere), so that the boundary of the BTZ string 
can be thought of as a sphere with a black hole at either pole. The bulk of the string is the 
interior of the sphere, where the string stretches from one black hole to the other. 



C-metric). This reference solution may be constructed by starting with global AdS4 writ- 
ten in coordinates for which slices of constant radial coordinate z are just AdSa. One then 
replaces each such AdSs slice with a BTZ metric \^ i33j having the correct 2;-dependence 
and which we chose to be nonrotating. The result is an AlAdS Einstein metric which may 
be written 



H^iz) 



-f{r)dt' + ^ + r'd<P' + dz' 
fin 



with H{z) = £3 cos^z/is) and /(r) = (r^ — rl)/i^. The solution is sketched in figure [2j Here 
the horizon of the BTZ string is at r = tq, the parameter £4 is the AdS4 length scale, and 
the AdSa length scale of the BTZ foliations may be set to any desired value by rescaling 
z,r,ro. The two boundary black holes (at z/i^ = ±7r/2) have the same temperature, but 
we will seek deformations where these temperatures differ and heat flows between the two 
boundary black holes. 

The outline of the paper is as follows. Section |2] reviews how static black funnels (i.e., 
without flow) may be generalized by adding a parameter a = Tbndy Bn/^buik bh which allows 
the temperature of bulk and boundary horizons to differ. In the small a limit, the analo- 
gous flowing funnels can be described in a derivative expansion; i.e., using the fluid/gravity 
correspondence of [31]. This correspondence is briefly reviewed and then applied to flowing 
funnels in section |3l Section |4] then explains how to formulate the construction of flowing 
funnels with any a in a manner where one can proceed numerically. The results of such 
numerics are presented in section O where they are compared to the fluid approximation of 
section [31 As expected, we find excellent agreement for small a, though for our cases the 
agreement remains good even for a close to 1. We close with some final discussion in section 

m 
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Note: In the final stages of this work we learned of |41j, which also addresses the con- 
struction of AdS black holes with non-Killing horizons and may have some overlap with our 
work. Their paper will appear simultaneously with ours on the arXiv. 



2 Detuning the bulk and boundary black hole temper- 
atures 

As mentioned in the introduction, even without heat fiow the black funnel paradigm may 
be generalized by adding an extra parameter a = J'"'^y which allows us to detune to 

-'bulk BH 

the temperatures of the bulk and boundary black holes. In terms of the dual field theory, 
taking a ^ 1 means that one considers a thermal ensemble at some temperature Tpieid Theory 
which differs from the natural temperature Tbndy bh of the (say, static) boundary black hole 
spacetime on which the field theory lives. One may think of the resulting state as defined 
by a Euclidean path integral with period l/Tpieid Theory 7^ 1/^bndy bh and thus having a 
conical singularity at the horizon of the boundary black hole. What is interesting about this 
construction is that the gravitational dual can have a completely smooth Euclidean AlAdS 
bulk, with the conical singularity of the boundary geometry resulting only from a failure 
of the standard AlAdS boundary conditions at the singular boundary points |l2l HSl 1^ . 
Any smooth horizon then clearly has temperature Tbuik bh = Tpieid Theory 7^ ^bndy bh as 
determined by the Euclidean period. 

The prototypical detuned solution studied in SI] is just the general hyperbolic 

(sometimes referred to as 'topological') black hole of SHlllZ], whose metric in D = d+1 
bulk spacetime dimensions may be written 

,2 T7l„\^-^2 , , 2j^2 TPI„\ ' 1 '0 / '0 



dsi,^, = -F{r)dt' + -— + r'd^U. = 7^ " 1 " (l^ " 1 ) • (2-1 



F(r) ^ ^ 



Here ^d+l is the AdS length scale associated with the [D = d + l)-dimensional cosmological 
constant, dTi\_^ = d^^ + sinh^ ^(if2^_2 is the metric on the unit Euclidean hyperboloid, and 
r = ro is a smooth horizon of temperature 

_ rlr,ld - (d - 2) 

-fbulk BH — • l-^-^J 

47rro 

Note that F{r) approaches r^/£^_,_j^ at large r. Making an obvious choice of boundary 
conformal frame, the boundary metric is just the hyperbolic cylinder H x M with (is^xR = 
—dt^ + C'^^j^^d'Ld-i- But note that we may write 

J 2 2 

where p/^d+i = tanh^. Multiplying the right-hand side by (1 — p'^/id+i) gives a metric on 
the static patch of the (i- dimensional de Sitter space dSd with Hubble constant id+i- ^ly 
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changing conformal frames in this way we may regard the boundary of (12.11) as having de 
Sitter horizons with temperature l/27r£d_|_i. From the perspective of an observer in the static 
patch, the de Sitter horizon acts just hke a black hole horizon with 



^bndy BH - —„ ■ (2.4) 

For general tq this temperature clearly differs from that of the bulk horizon. For the case 
where they agree, the hyperbohc black hole metric (12. ip is just pure AdS^+i in appropriate 
hyperbolic coordinates. We recall that even for the tuned case a = 1 ref. [21] found the 
conformal frame (12. 3p useful for constructing black funnel solutions numerically. 

Since the analysis of temperatures above involves only the horizons, it is clear that 
detuned bulk and boundary horizons should exist much more generally. Indeed, any static, 
spherically symmetric boundary metric with a pair of of smooth horizons at p = ±£^+1 may 
be written in the form 

dsU, BH = (1 - pV^^+i) [--nP)dt' + 1^ + R\p)d^l-^ ' (2-5) 

where F, G, and IjR^ are smooth on some interval including p G [—(■d+iAd+i\, G has a 
second order zero at each of p = ±£^+1, and 1/B? vanishes at p = ^dd+i- So after a 
conformal transformation (12. 5p agrees with (12. 3p to leading order in p for each term and in 
this sense may be said to approach H x M at large p. The ansatz (12. 5p can equivalently be 
written 

^^^Ldy BH = e-2-/-oF(s) {-de + dx^ + R\x)dnl_^) , (2.6) 

where Xq is some reference length scale and F and e^^^/"^"/?^ are smooth functions of e^^^/^o 
at e^^^/^° = 0. In particular, up to the conformal factor rfs^xR takes this form for xq = id+i 
and = sinh^(a;/£d+i)- In terms of (12. 6 p the boundary black holes have temperatures 

ybndy BH ^ ^ d / -2x/xo^/^)N _ ^2.7) 

47r x^±oo dx ^ ' 

It is therefore sensible to choose any ro and seek a smooth bulk solution in which each 
term approaches that of (12. ip to leading order in at large see section H] for a 

more complete analysis of these boundary conditions. Any static such solution will have a 
bulk horizon with temperature (12. 2 p and can again be interpreted as being dual to a field 
theory state of this temperature on a black hole background of temperature (12. 4p . In the 
next sections we will seek further generalizations with different values of ro (which we then 
call r±) at x = ±oo. That is to say that for x — )■ +oo the bulk solution will asymptote 
as above to (12. ip with ro = r+, while for x — )■ — oo it will analogously approach (12. ip with 
ro = r_. The bulk horizon may then be said to approach the temperatures T± given by 
(12. 2p with ro replaced r±. We will also allow distinct temperatures T^^*^^ for the x = ±oo 
boundary black holes and introduce the parameters a± = T^^'^^ ^^/T±. In fact, we will 
always take a+ = a_. 
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Of course, we may also consider so-called ultrastatic conformal frames analogous to (12.31) . 
Starting with (12.61) and multiplying by a conformal factor e^^'/^'°/F(a;), one obtains the 
boundary metric 

ds'^ = -dt^ + dx^ + R^{x)dnd-2 (2.8) 

for which dt is a hypersurface-orthogonal Killing field of norm —1. In this frame, the bound- 
ary spacetime has two asymptotic regions, each asymptotic to H x R (say, with the same 
curvature scale £d+i)- Furthermore, in the CFT description each region contains an infinite 
reservoir of deconfined plasma. Such infinite reservoirs may act as heat baths, and indeed 
the boundary conditions imply that they are in thermal equilibrium at temperature T± in 
the limits x — ?■ ±oo. 



3 The fluid limit 

While a general treatment of black funnels remains challenging, it is by now well known 
that the study of AdS black holes simplifies in the so-called hydrodynamic limit of the 
fluid/gravity correspondence j34j in which all other parameters vary slowly in comparison 
with the black hole temperature and the solution can be described using a derivative expan- 
sion. For any fixed boundary metric, taking the limit of large temperature (i.e., small a±) 
makes all metric derivatives small in this sense. We may thus expect a good hydrodynamic 
description if in addition we control temperature gradients by taking AT = T+ — T_ small. 

The key point in the analysis of [34j is that, having chosen a boundary conformal frame 
with boundary metric glj\ every AlAdS^+i solution is associated with a (i-dimensional 
boundary stress tensor Tij which is traceless and conserved on the boundary: 

glfT^ = 0, Dir^ = 0, (3.1) 
where Di is the covariant derivative compatible with g^j"^ . Below, we use the boundary 



metric gf^^ and its inverse to raise and lower indices k,l, . . . . 



As an example, consider the planar AdS-Schwarzschild black hole 
dsl,s-s... = - {rViU (1 - de + 1^%,^^,^ + {rViU d^U, (3.2) 

with To = 4:7ii'^_^_j^T/d. Taking the boundary metric to be ds^^^y = —dt^ + (ix^_^, one finds 

T''=TiL^ = P^^+PV'^. (3.3) 

which takes the form of an ideal fluid with velocity field u^di = dt, transverse projector 
= gij + u'u^^ and 

which of course satisfies fl3.ip . In (13.41) . we have defined for convenience T = ini^^iT/d. By 
a simple Lorentz transformation we may obtain corresponding solutions with any constant 
(normalized) timelike rf- velocity u\ 
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The main result of |3l] was to show that the temperature T and d-velocity may be 
promoted to slowly- varying functions of the boundary coordinates x, t (at which point we 
refer to them collectively as the hydrodynamic fields). Here the term "slowly- varying" is 
defined with respect to the temperature as measured in the local rest frame selected by 
ti*. In particular, under these conditions [M] showed that a smooth bulk solution may be 
constructed via a gradient expansion so long as is everywhere timelike and the associated 
boundary stress tensor does indeed satisfy fl3.ip . They further showed that at each order in 
this expansion the conditions fl3.1l) may be expressed as standard hydrodynamic equations 
for a (conformal) fluid with velocity field m*, which we take to satisfy u^Ui = —1. This last 
step essentially just repeats the standard derivation of hydrodynamics from conservation 
laws. 

In particular, ref. [51] showed that the boundary stress tensor takes the form 

T''=T1^ + Y.^^ (3-5) 

n=l 

where are dissipative terms that are n**^ order in derivatives of the hydrodynamic fields; 
for example, 

n;^,) = -2va'\ (3.6) 

where 

is the shear viscosity, and 6 = DiU^ and 

^ii = pi'^VJiD^kUi) - ^ r^' (3.8) 

are respectively the divergence and shear of the velocity field. In writing fl3.6p there is a 
freedom to make certain field redefinitions which, following [51], we have removed by choosing 
the so-called Landau frame in which the n*^-^-| are taken to be purely transverse. 

Since by assumption derivatives of the hydrodynamic fields are parametrically small in 
some parameter e, !!(„) is of order e". Below, we solve the fluid equations (13. ip at order 
n = and n = 1 for the ultrastatic boundary metrics (12. 8 p and a purely radial velocity 
field (so that the only non- vanishing components are m*, u^). We also assume the flow to be 
stationary, so that u*, T are independent of time. 

A new effect at first order is the appearance of dissipation, and thus the production 
of entropy. At zeroth order, the entropy current Jg takes the simple form {Jg)idea.i = su\ 
where s{x) = T^~^ / AG^"^^^^ is the entropy density. Using the equations of motion and 
thermodynamic relations, one can show [48j that 

AW),aeai = 0- (3-9) 
At first order, the entropy current still takes the form (^5)1 = su^, but its divergence now 
becomes [H] 

A(4), = ^^a,a«^>0, (3.10) 
showing that entropy is produced unless aij = 0. 
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3.1 Ideal Fluid 



We begin at order n = 0. We denote the associated fluid quantities %, Mq and work in c? = 3. 
Following [19], we project the fluid equations into components parallel and perpendicular to 
the velocity. These yield respectively 



or 



Thus 



A (To^'o) = and D'To + n^A {To4) = 
■ and {% {uo)^] 



d. 



0. 



•2 

oo 



2aR' 



To (Mo)t = Toe 



(3.11) 
(3.12) 

(3.13) 



in terms of integration constants that we have chosen to call T^/2a, Too- Since = —1, it 
remains to solve a quadratic equation for To,Uq. We of course obtain two solutions labeled 
by a choice of sign. The solution with finite and nonzero asymptotic temperatures T± has 



~ 2 



1 - 



Un 



aR 



(3.14) 
(3.15) 



Note that since R diverges at large x, at this order the asymptotic temperatures T± at 
X — 7- ±oo agree; i.e., AT = — = 0{e). We also find — > at a; = ±oo. 



3.2 First Order Corrections 

To compute corrections to fl3.14p . fl3.15p . we choose to solve the fluid equations fl3.ip iter- 
atively. Introducing a bookkeeping parameter e to keep track of derivatives, we may write 
T = Tm + 0{e^^^), = + C(e™'+^) for each m. We compute 7^, by dropping terms 
with n > m in (13. 5 p and evaluating the remaining n^-^^ on Tm-n, ^m-n- Thus Tm, u^m enter 

(13. ip only through T^^eai equations to be solved are essentially just (13. lip . (I3.12p 

with additional source terms given by the H*^-^^. The integration constants (as well as the 
sign choices that come from solving quadratic equations) may be flxed by requiring 7^, 
to approximate Tm-i, ^m-i ^^e desired order as e — )■ 0. 
To flrst order, one flnds 



T^ = \{B{x)+Toof 



{B{x) + ToofR 
2(A(a;) + ri/2a) 



1 + 



( 2{A{x) + rma) 
V {B{x) + TooYR 



(3.16) 



( 2{A{x) + rma) 
V {B{x) + TooYR 



(3.17) 
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where 



3 Jo 



B(x) 



24 
~3~ 



(0) 



(3.18) 
(3.19) 



and the square brackets [-j*^"^ indicate that the enclosed quantities are evaluated on the zeroth 
order solutions fl3.14p . f l3.15p . At this order, the asymptotic temperatures differ and are given 
by the (finite) expression 

r(±oo) = roo + 5(±oo), (3.20) 

so that 

AT := r(oo) - r(-cx)) = 5(00) - B{-oo). (3.21) 

It is useful to consider the further limit of small AT, which greatly simplifies the above 
results. This is equivalent to taking a large. Since 

2" 



Bix] 



2^4 



R"{x') 1 



2/22 (a;') a 



+ 



dx'. 



we find 



for 



2R'^(x') a 



dx' = -— + O i - 



R"{x) 
RHx) 



dx. 



Noting that A{x) = 0{l/af we then find 

u\ = i + o{Ar^), 

SAT 



Ml 



Ti = Too + 



2IR{x) 

\Ar 



^ R"ix' 







^ dx' + o(Ar' 



so that the non-zero components of the stress tensor are 

QhJlAT r R"(x'^ 



-2T: 



•3 



R^ix 



^ dx' + 0{AT' 



QT^AT 



2IR{x) 



UaJIAT fR' 



[x] 



RHx) 



R'^ix' 
" R"(x' 



R^ix' 



jrdx']+ 0{AT' 



(3.22) 

(3.23) 
(3.24) 

(3.25) 
(3.26) 

(3.27) 

(3.28) 
(3.29) 
(3.30) 
(3.31) 
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Note that the lowest order term in the energy flux T^^ is linear in AT; this naturally 
leads to a notion of thermal conductivity. We first calculate the heat flux $ as the energy 
flux integrated over a circle of constant x: 



$ = 2txR{x)T 



tx 



QT'^ AT 

— + c>rA'D 



(3.32) 



We define the thermal conductivity as := — (i$/(iAT|AT=o so that 

k ^""^^ 



4^(4) J ■ 



(3.33) 



We have also explored the analogous results at second order n = 2 m the hydrodynamic 
approximation. While the general expressions are unenlightening, each quantity above agrees 
with the n = 1 expression up to linear order in AT for all Too- In particular, the conductivity 
k is unchanged. 

Finally, the entropy current {Js)i = sm* for our solutions is 



{B{x) + rooYR 

2V2{A{x) + ri/2a) 



2{A{x) + ri/2a) 
{B{x) + TooYR 



1/2 



X 1 + Wl 



2{A{x) + Tl/2a) \ 
{B{x) + TooYR J 



(3.34) 

(3.35) 
(3.36) 



which has divergence 



R''- 



R^{a^R^ - 1) 

To lowest nonvanishing order in AT, these become 

2£4TooAT 



a?R^ 



1/2 



(3.37) 



4G(^)(4)^=T^+ ^ 

9^-2 A T 

4G<'" ( J|), = 



2IR{x) 

AG^^^D, (4)^ = MlZ^ + o(AT)= 



(3.38) 
(3.39) 

(3.40) 



P 

Note that the divergence of the current is of order AT^ as expected from fl3.10p . It turns out 
that fl3.40p is unchanged when one passes to second order in the hydrodynamic expansion, 
though the entropy current Jg itself changes even at zeroth order in AT. 
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These expressions may of course be transformed to any other conformal frame. The 
ultrastatic frame fl2.8p used above had the convenient feature that, at least at small velocity, 
the local fluid temperature (defined with respect to proper time in the fluid rest frame) 
coincides with the temperature defined with respect to the static Killing field dt. In a more 
general conformal frame, these two temperatures do not coincide even at small velocity. 
Note that we will employ only time-independent conformal transformations below, so that 
dt remains a Killing field in all frames. We will continue to refer to temperatures normalized 
(up to a boost to the fluid rest frame) with respect to dt by T, while we denote the local fluid 
temperature (defined with respect to rest-frame proper time) as TToc- Thus T is unchanged 
by the conformal transformation while TToc is rescaled. 

For comparison with our later numerics, appendix |X] presents the results in the black hole 
frame for the explicit metric functions and in terms of the particular coordinates used in sec- 
tion E] below. The resulting more explicit expressions are correspondingly more complicated 
than those above. 



4 How to flow a more general funnel 

Our family of flowing funnels will be labeled by four parameters: the temperatures T^"^^ 
of the left- and right- boundary black holes and the temperatures T± associated with the left- 
and right- ends of the bulk black hole. As discussed in section H] these four temperatures are 
completely independent in principle, though in our simulations we will always set q;+ = a_ 
which introduces one relation. 

The most generic ansatz compatible with our symmetry requirements depends on seven 
unknown functions: 



1 — w^Yil — y 



2\2 



-M{y)G{w)\l-wyy' A 



?J^di+Q(w)—dy 

y 



+- — TTT^ — - + yl 



4 5i 



dw + £4 ^xi dt + 



Fdy 



+ S2 



, (4.1) 



M{y) '''' \ - " y 

where A, B, F, Si, S2, Xi and X2 are all functions of y and w. In addition we have defined 

2 



G(w) = l+-w\5-3w^) , M(y) = 2-y^-^ — ^^^^^ — ^ and Qiw) = 1+ 
2 2/0 



M{0)G{w) 
(4.2) 

The insertion of these factors will be justified later, when we will also see that /3 controls 
the temperature difference between the two boundary black holes, and yo is a parameter 
that controls the validity of the fluid approximation. Here y ranges over [0, 1] and w ranges 
over [—1, 1], with y = being the bulk horizon and y = 1 the conformal boundary. At least 
for y ^ regions with w ~ ±1 are close (in the sense of a conformal diagram) to where 
either bulk horizon meets either the left or right boundary black hole (compare with figure 
[2]). As we will explain below, the symbol t was used in (14. ip in order to save t for another 
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coordinate associated with Feffer man- Graham gauge. However, di = dt so we will refer to 
the time-translation as simply dt- 



4.1 Boundary Conditions 

At the conformal boundary (?/ = 1) we impose the boundary conditions 

A{w, 1) = B{w, 1) = Si{w, 1) = S2{w, 1) = X2{w, 1) = 1, F{w, 1) = xi{w, 1) = 0, (4.3) 
which ensures a boundary metric conformal to 

ifdsl = -^(1 - pYGipr dt^ + ^ + #2 , (4.4) 

where p = p/i4^. As in section |2l we refer to (14. 4p as the boundary metric in the black 
hole conformal frame. In presenting our results in section E] we will describe all boundary 
quantities, such as the stress energy tensor, with respect to this frame. The boundary metric 
dsg has horizons at p = ±1 with Hawking temperatures 

^bndy BH _ G{±1) 

We will extract the boundary stress tensor following the strategy of [31] and using the results 
of [50]. The only technical difference with respect to [31] involves the relation between the 
coordinates (t, w, y, 0) and Fefferman-Graham coordinates {t, z, p, 0). Due to the cross term 
X2 in Eq. (14. ip the map between i and t is not trivial, instead it is expressed as a powers 
series in z of the form: 

i = t + zTi{p) + 0{z^) (4.6) 

where for instance Ti(p) = —(^(p) 2/0/ (2(1 — p^)). 

The left and right boundaries lie at = ±1. There we impose 

A{±l,y) = B{±l,y) = Si{±l,y) = S2{±l,y) = X2{±l,y) = 1, F{±l,y) = Xi{±hy) = 0, 

(4.7) 

which reduces Eq. (14. ip to 



?rMt + g(±i)- 
y 




Under the coordinate transformation: 



y=f^^ |c^r = G(±l)(£-cit+^i^), rv = ±lTe 



1/0 = ?, (4.9) 
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the line element f l4.8p yields the large ^ limit of Eq. (12 .ip with d = 3. The fact that our 
ansatz fl4.ip reduces to a hyperbolic black hole at w = ±1 displays the physical meaning of 
yo as an overall scale that controls the bulk horizon temperatures (and thus a±). Note that 
the line element (gSD also defines T± = T^"''^^ ^"M(0)/2. If yo = 1, then T± = T^"'^^ i.e. 
it represents the 'tuned' case a± = 1. Thus the fluid approximation becomes more accurate 
as yo increases, or equivalently, as a± decrease. 

We have imposed Dirichlet data at each of the above three edges of our computational 
domain. But it remains to specify boundary conditions at y = 0, the flowing funnel horizon. 
Here we demand that the line element (14. ip be smooth in ingoing Eddington-Finkelstein co- 
ordinates (which cover the future horizon). To understand the explicit form of this condition, 
we introduce local ingoing Eddington-Finkelstein coordinates {v,w,y,(f)) through 

dv = dt + U% + 0{y^), dw= f"^^. +i^'dv + Oiy''), y = y^'\ (4.10) 

The terms omitted in the above y expansion can be chosen such that a line of constant 
{v,w,(j)) is an ingoing null geodesic. Note that lines of constant v have dy/di < 0, as 
required for ingoing coordinates. Furthermore, regularity of the line element (14. ip in the 
above coordinates requires 

F(., 0) ^ xi(^, 0) , Bi., 0) ^ ^(^W^(-^0) ^^^^^^^^^ 

dyA{w,0)=0, dySi{w,0) = 0, dyS2{w,0) = 0, dyXi{w,0)=0, dyX2{w,0)=0. 

(4.11) 

We will find Xi(^) 0) to be finite and non-zero (at w ^ ±1), so our original w is already an 
ingoing coordinate. It will thus be straightforward to read off results associated with the 
future horizon. 

The past horizon is more subtle. It is located at v ^ — oo and can be reached along lines 
of constant w. Depending on the sign of xi, this tends to drive w to either ±1. Below, we 
consider > T_ so that the hotter black hole is on the right. One might therefore expect 
w to decrease along the horizon generators so their coordinate velocity is toward the cooler 
black hole; i.e., one might expect Xii'^y 0) > 0. But for the particular ansatz we have chosen 
our numerics turn out to give Xi < (see section E]) so that the past horizon in fact lies at 
w = —1. This appears to be a coordinate artifact, though a full understanding is beyond 
the scope of this work. 

Below, we will solve the Einstein equations (with cosmological constant) in the form 

Eab := Rab + igab = 0, (4.12) 

subject to the boundary conditions detailed above. 



4.2 The DeTurck Method 

The diffeomorphism invariance of fl4.12p means that these equations do not lead to a well- 
posed boundary value problem. While one could attempt to proceed by gauge-fixing, a 
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clever trick known as the DeTurck method was introduced in [27] and in [28j was shown to 
succeed (under rather general assumptions) when one seeks appropriate stationary equilib- 
rium solutions of the vacuum Einstein equations, with or without a negative cosmological 
constant. Though our situation turns out to fall outside the bounds of the proof given in 
[28] . we nevertheless employ this method successfully below. 

We begin with a brief review. The DeTurck method is based on the so called Einstein- 
DeTurck equation 

E^, = - V(,4) = 0, (4.13) 

which differs from from Eq. (14.121) by the addition of —V(atb)- Here = g'^'^\r^^^{g) — T'^^^{g)], 
is the Levi-Civita connection associated with the metric g, and g is some specified non- 
dynamical reference metric. Since ^ is defined by a difference between two connections, it 
transforms as a tensor. Hereafter g will be chosen to have the same asymptotics and horizon 
structure as g. In particular, it must satisfy the same Dirichlet boundary conditions as g. 

Clearly any solution to E^^ = with ^ = also solves Eat = 0. But one may ask if (I4.13P 
can have additional solutions that do not satisfy Eab = 0. Under a variety of circumstances 
one can show that solutions with ^ 7^ 0, the so called Ricci solitons, cannot exist [25] . 
However, the assumptions used in [2S] seem not to hold for our system of equations. In 
particular, after reduction along the symmetry directions t, our system turns out to have 
a mixed-elliptic hyperbolic nature. This is most easily seen from the fact that, while our 
system will be elliptic near infinity where dt is timelike, we expect an ergoregion near the 
horizon where all linear combinations of dt, are spacelike. So in this region reduction 
along (t, 0) gives a Lorentz-signature metric on the base space. This differs qualitatively 
from the case of Kerr, where dt-, span a timelike plane everywhere outside the horizon and 
reduction along (t, 0) gives a Euclidean-signature metric on the base space. The difference 
arises from the fact that the Kerr horizon 'flows' only along the Killing field while our 
horizon 'flows' in the w direction, which is not associated with any symmetry. Thus Ricci 
solitons may well exist in our case. But for any solution to (I4.13P one may simply calculate 
to see if it vanishes. For all of our flowing funnel solutions discussed below we find ^ = 
to machine precision. 

It remains to specify our choice of reference metric g. We choose g to be given by the line 
element (14. ip with A = B = Si = S2 = X2 = ^ and F = xi = 0. This enforces all Dirichlet 
boundary conditions except those at the horizon, Eq. (14. lip . To satisfy these remaining 
conditions we need only choose Q{x) as in Eq. (14. 2p . 

4.3 Numerical Method 

We use a standard pseudospectral collocation approximation 'm.w,y and solve the resulting 
non-linear algebraic equations using a damped Newton method with damping monitoring 
function This ensures that Newton's method takes a path in the approximate solution 
space that decreases \S,t\ at each step. This method may also prove useful in solving more 
general mixed elliptic-hyperbolic systems. We represent the w and y dependence of all 
functions as a series in Chebyshev polynomials. As explained above, our integration domain 
lives on a rectangular grid, {w,y) G [—1, 1] x [0, 1]. 
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Figure 3: Ajy as a function of the number of grid points N. The vertical scale is logarithmic, 
and the data is well fit by an exponential decay: log(Ajv) = —17.4 — 0.23 A^. 



To monitor the convergence of our method we have computed the total heat flux $ 
(defined by the first equality in (13.321) ) for several resolutions. We denote the number of grid 
points in w and y hj N and compute Ajv = |1 — ^n/^n+i\ for several values of A^. The 
results for this procedure are illustrated in Fig.Ofor /3 = 0.1 and = 1. We find exponential 
convergence with A^, as expected for pseudospectral collocation methods. Furthermore, in 
order to ensure that we are converging to an Einstein solution rather than a Ricci soliton we 
monitor all components of ^. For all plots shown in this manuscript, each component of 
has absolute value smaller than 10^^°. 



5 Results and comparisons 

We now present the results of our numerical analysis and compare them with the first-order 
(n = 1) hydrodynamic approximation. The plots below are labeled by a parameter Too 
whose definition 

_ [256 (144^2 - 557) + 1057r (293A2 + 128)] (3?/g - 1) 
~ ~ 287r [157r (293A2 + 128) - llOOSA^] y^ ^ ' ' 

was inspired by the first-order hydrodynamic result (I3.16p . For small AT we have Too = 
(T+ + T_)/2 + (9(AT)^. We note that all the numerical results we will present use units 
where £4 = 1 (so that p = p) and 16nG^^^ = 1. We also take T+ > T_ so that the hotter 
black hole lies on the right. 



We begin with the norm |(9tp of the time translation. Figure 4(a) shows a typical plot. 
To guide the eye we have also plotted a reference surface of constant |(9tp = 0. The two 
surfaces intersect at the ergosurface, whose location we display separately in Fig. 4(b) Inside 
the ergoregion |c}tp becomes positive, changing the character of Eq. (I4.13P from elliptic to 
hyperbolic. This region is at the core of the difficulties in trying to prove that our numerical 
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(a) (b) 

Figure 4: (a): The curved surface shows the norm of dt over our integration domain. To 
guide the eye, we also draw a flat horizontal surface at zero norm, (b): The ergosurface as 
a function of w. Both figures use a± = 1 and AT/Tqo = 0.2. 

method ensures ^ = on solutions of Eq. f l4.13p with appropriate boundary conditions. Fig. |5] 
shows \dt\'^ and, for comparison and later use, |9<^p as a function of w along the horizon. We 
remind the reader that dt and are precisely orthogonal everywhere in our spacetime, so 
this describes the full induced metric hjj (for I, J = t, 0) in the 2-plane spanned by dt, d^. 
Both norms are clearly positive everywhere on the horizon, though \dt\'^ never becomes very 
large even with AT/Tqo = 0.2. This may help to explain why our numerical approach 
succeeded. 

Let us now discuss the behavior of the boundary stress tensor. For small a, AT/T^, 
this quantity may also be computed using the hydrodynamic approximation of section [3] and 
provides another good check of our numerics. Fig. [6] shows the components of the stress 
energy tensor as a function of the boundary coordinate p for several values of a at fixed 
p. The lines represent the first order hydrodynamic prediction and the symbols represent 
data extracted from our numerics. Large stress tensors correspond to larger values of a. 
We see that at least for small AT/T^ the fluid gravity prediction works remarkably well 
even for for a ~ 0.8. The agreement of all of these curves when a is small is a reassuring 
test of our numerics. However, at larger a qualitative differences from our hydrodynamic 
approximation begin to appear. For example, we note that while T*f is always negative (and 
thus the energy density is positive) in the hydrodynamic regime, for a > 1 our simulations 
show T*t becomes positive near the hotter black hole. 

From the standpoint of the dual CFT, the main physical result of our paper is displayed 
in Fig. [71 This plot shows how the total heat flux $ varies for different values of AT /Too 
and a = a+ = a-. We see that it increases in magnitude as AT/T^o increases, and also as 
a decreases. This computation can be seen as a first principle calculation for the thermal 
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Figure 5: (a): The norm \dt\'^ on the future horizon, (b): The norm \d(j,\'^ on the future 
horizon. Both figures use a± = 1 and AT/Too = 0.2 and are plotted as functions of w. 



conductivity of a strongly coupled plasma at large beyond the hydrodynamic regime. 
Fig. |8] compares some a = const, cross-sections of Fig. [7] to the the results of first-order (n = 
1) hydrodynamics at linear order in AT; i.e., to flA.2p -f lX3|) . These show good agreement 
for small a and AT, but deviate as expected at larger a. 

It remains to examine the horizon more closely. Our horizon is a three-dimensional 
null surface and, since dt, d,p are both spacelike and tangent to the horizon, any two null 
geodesies that generate the future horizon generators are related by some isometry. Thus all 
generators are equivalent, though it remains to understand the evolution of the spacetime 
along each generator. We compute the affine parameter, expansion, and shear along each 
generator using simple expressions in terms of the induced metric hjj (for I, J = t, 0) on 
the 2-plane spanned by dt,dij,. These expressions are given in appendix [Bl We study each of 
these quantities only on the surface y = 0. 

We begin with hjj itself. Recall that w = ±1 are the asymptotic regions of static 
hyperbolic black holes where the Killing field dt becomes null at the horizon and |9<^p becomes 



large. These behaviors are clearly shown in figure 9(a) But these similarities between 



w = ±1 are misleading and the actual behaviors at w = ±1 are quite different. This may 



be seen from the plot oi h = det hjj = huh^f,^ = \dt\ \dfj,\ in Fig. 9(a) This determinant 
vanishes at w = — 1 but approaches a non-zero constant aXw = +1. Note that h is monotonic 
along ?/ = 0, as it must be along a smooth horizon. 

What is perhaps surprising is that h is an increasing function of w. This shows that 
w increases toward the future along the future horizon, so that the past horizon must lie 
at w = — 1. In contrast, in the coordinates of e.g. |51J, the coordinate velocity of the 
horizon generators would be in the direction of heat transport, and thus (since we take the 
cooler black hole to lie a.t w = —1) toward negative w. Standard coordinates for Kerr also 
behave like those of ^T] and have the equivalent of our xi being positive for positive angular 
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Figure 6: Components of the boundary stress energy tensor as a function of p for fixed 
/3 = 0.04. Each panel shows the first-order {n = 1) hydrodynamic prediction as hues 
and the exact numerical data as symbols. The disks and solid line show a-t = 1, the 
squares and dashed line show a± = 0.77 and the diamonds and dotted line show a± = 0.70. 
These corresponds to /ST /Too = 0.080, AT/Tqo = 0.050 and AT/Tqo = 0.034, respectively. 
Since AT /Too is small, we have used only the linear results from appendix |X] to plot the 
hydrodynamics. 



velocity. In contrast, we find xi to be negative at the horizon; see figure 9(b) Since xi 



samples completely different metric components than h, we take this as a strong indication 
that our solutions are consistent despite the surprising location of the past horizon. Another 
strong indication of consistency is the above agreement between our boundary stress tensors 
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Figure 7: Three-dimensional plot of the boundary flux extracted from our numerics as a 
function of AT /Too and « = = 




" •] |_ 
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Figure 8: The total heat flux <l> as a function of AT/T^o for a = 0.9 (left) and a = 0.7 
(right). The solid curves are the first order hydrodynamic results. Since AT /Too is small, 
we have used only the linear results from appendix |Al The dots show our numerical data. 

and those predicted by the hydrodynamic approximation. Indeed, we have tested for various 
possible errors (such as inverting the sign of AT) in our code by examining the effect of 
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Figure 9: (a): The determinant h = huh^^. (b): The metric component Xi along the horizon. 
Both quantities are plotted as functions of w for a± = 1 and AT /T^o = 0.2. 



various sign changes on Fig. [H] and found in each case that such changes would lead to 
notable discrepancies with hydrodynamics. In particular, we stress that our simulations give 
the physically correct sign for the heat flux T' 



The apparent proximity of the past horizon to the cooler black hole must thus be a 
coordinate artifact. While we would like to understand this effect in more detail (say, by 
changing to the coordinates of [5T] for a solution deep within the hydrodynamic regime), 
the additional numerics required places such an analysis is beyond the scope of this work. 
We note, however, that a similarly surprising sign can be found in the 2-1-1 flowing funnels 
of [30]. In that case, writing the horizon generating Killing field in the Fefferman-Graham 
coordinates of [30] leads to a negative t component on part of the horizon, even though this 
component is everywhere positive at the AlAdS boundary. 

We may now proceed to investigate various quantities along the horizon. Perhaps the 



most important quantity is the affine parameter A, which we show in Fig. 10(a) as function 
of w. Note that A approaches a constant value at w = —1. This is to be expected, as we have 
already noted that w = — 1 is the past horizon. Since the affine parameter is only defined 
up to affine transformations, this constant is arbitrary and we have set A(tf = — 1) = for 
convenience. In contrast, the affine parameter diverges as we approach w = 1. 



Figure 10(b) shows the expansion 6^ as a function of A. As expected on general grounds. 



9 is everywhere positive with d6/dX < and 9 asymptotes to zero at large A. We see this as 
the most solid test of our numerics. Note that the sign of d9/dX < is only guaranteed via 
Raychaudhuri's equation once the equations of motion are used. It is thus far from trivial 
that the sign comes out right. 

The expansion diverges at the past horizon (A = 0), indicating the presence of a caustic. 



In fact, it is easy to see that this caustic is a curvature singularity. To do so, note from 5(b) 
that \d(t,\ diverges on the past horizon. But since Killing fields obey a second order differential 
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Figure 10: (a): An affine parameter along the horizon as a function of w. (b): The expansion 
of a future horizon generator as a function of A. (c): The positive eigenvalue a of 
function of A. At small A we find a ~ X~^^^. All figures use a± = 1 and AT/T^c = 0.2. 



equation governed by the Riemann tensor (see e.g. (C.3.6) of fS2]) they can diverge at finite 
affine parameter only if Rated diverges in all orthonormal frames. 

We now turn to the shear tensor From ( 1B.4|) . ( 1B.7|) we see that since hu is diagonal, 
the same is true of a/j . Since ajj is also symmetric and traceless, it is completely charac- 
terized by the positive eigenvalue cr of a^j, where the index was raised with the inverse 
of hjj. Note that a is a spacetime scalar. 

As one might expect, it 



This eigenvalue is plotted as a function of A in figure 10 (c 



diverges at the caustic. What is interesting is that we find the same divergence structure for 
all a, AT that we have studied. We quantify this behavior by fitting o"(A) to a power law 
IJ,X^ near A = 0. We have extracted rj for about 400 fiowing funnels spanning the domain 
(a, at/Too) e (1,0.7) X (0,0.4). In all cases we find t] = -0.82 ± 0.03, where this error is 
in fact the maximum deviation. We note that this number is remarkably close to —5/6. We 
can then use the Raychaudhuri equation f IB.Sp and the standard evolution equation for the 
shear (see (F.34) of |53]) to again show that Rabcd diverges on the past horizon. In fact, for 
7] = —5/6 one may show that some Weyl tensor component Cabcdk^k'^ (where fc" is an affinely 
parametrized tangent to a null generator of the horizon) must diverge like A~^^/^. This fact 
merits an analytic explanation which we are unable to offer at this time. 



6 Discussion 

The above work constructs 'flowing funnel' stationary black hole solutions. Such solutions 
describe heat flow between reservoirs at unequal temperatures T±. The particular solutions 
constructed are global AdS4 flowing funnels which may be thought of as deformations of the 
BTZ black string (11. ip . Thus each heat reservoir lies just outside a boundary black hole 
of temperature a±T±. For the case a± = 1, the CFT3 duals of our bulk solutions describe 
heat transfer between two non-dynamical 3-dimensional black holes due to CFT3 Hawking 
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radiation . 

Our solutions display many properties expected on general grounds. There is a connected 
ergoregion near the horizon where dt becomes spacelike. In fact, all Killing fields are spacelike 
at the future horizon "H, so that "H is not a Killing horizon. This is consistent with the rigidity 
theorems [371 EH El] since "H is not compact. 

The expansion 9 of the null generators is everywhere positive but decreases toward the 
future along each null generator. The generators extend to infinite affine parameter in the far 
future (where 6* — )■ 0) but reach a caustic {6 — )■ oo) at finite affine parameter toward the past. 
This caustic occurs on the past horizon, which is a curvature singularity characterized by a 
universal power law divergences for the shear a ~ A~^/^ and for certain components of the 
Weyl tensor which grow like A~^^/^ in any orthonormal frame. These exponents were found 
numerically, but merit an analytic understanding. It remains an open question whether 
curvature scalars (e.g. the Kretschmann scalar Rahcd.R"'^'^'^) niight remain finit(EI. It would 
also be interesting to study the exponents governing the divergence of the expansion 9, the 
norm jc^^p and the inverse norm though these have proved to be more difficult to 

extract from our numerics. 

Note that \d^\ decreases with A along the early part of the future horizon. But since 
^ > 0, the shrinkage of the circle with affine parameter A is more than compensated by 
the growth in \dt\. This positive expansion is associated with the expected generation of 
entropy due to the transport of heat from a hot source to a cold sink. In particular, it is the 
analogue at large a±, AT /Too of the entropy generation term fl3.10p seen in the hydrodynamic 
approximation. 

In our coordinate system, the horizon generators appear to fiow toward the hotter black 
hole. As our results pass a number of consistency checks, this seems to be a coordinate 
artifact. It would nevertheless be desirable to understand this effect in more detail. 

We studied the boundary stress tensor of such solutions both numerically and to first 
order in the hydrodynamic (fiuid/gravity) approximation. In particular, we computed the 
total heat fiux $ for boundary metrics of the form f l4.4p as a function of a, AT/T^o] see figure 
[71 It would clearly be of interest to study more general boundary metrics to understand which 
parts of this function are universal and which depend on detailed features of the boundary 
metric. 

The hydrodynamic approximation is a derivative expansion which, since we fix all other 
parameters in the boundary metric, is for us governed by the parameters a± (which con- 
trol the extent to which the bulk and boundary black hole temperatures are detuned) and 
AT/ (T+ + T_) (which controls the temperature difference between the heat source and sink). 
As expected, we find excellent quantitative agreement when these parameters are small. In- 
terestingly, we also find good qualitative agreement when these parameters are close to 1. 
This gives yet another confirmation of the robust nature of the fiuid/gravity correspondence 
as seen previously in e.g. [3j. Of course, at large enough values of a±, AT/ (T+ + T_) we find 
both quantitative and quahtative disagreement. It would be interesting see to what extent 
agreement might be improved by incorporating higher order hydrodynamic corrections. A 

^Due to the large number of terms involved, we were not able to reliably calculate RabcdR"'^^'^ from our 
numerics. 
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particularly notable feature at large a (a+ ^ 1 in our simulations) is that, while T*^ is always 
negative in the hydro dynamic limit, it becomes positive close to the hotter boundary black 
hole black hole. It would be interesting to understand this feature analytically. 
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A Fluid results in the black hole frame 



We may transform the hydrodynamic results of section 13.21 to the black hole frame associ- 
ated with the metric f l4.4p by implementing a boundary conformal transformation and an 
appropriate change of coordinates. Setting £4 = 1 the result is 
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where AT is again defined with respect to dt, Tioc is the local temperature with respect to 
proper time in the fluid rest frame, and setting H{p) := (1 — p^)G{p)/yQ we have defined 



h{p) = ^-./2^H{p)H'{p), 
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B The horizon-generating null congruence 



We wish to study the expansion and the shear tensor associated with the null geodesic 
congruence that generates the future bulk horizon. Instead of solving the geodesic equation 
and taking derivatives of deviation vectors, we take advantage of the fact that our system 
is co-homogeneity 2 to compute these quantities directly (up to a position-dependent scale 
factor) from the induced metric hjj on 2-dimensional surfaces tangent to the Killing fields 
dt, d^. We then compute the affine parameter A along these geodesies from the Raychaudhuri 
equation as explained below. 

Recall that we consider a future event horizon "H of an AlAdS4 spacetime with two 
commuting KVFs dt and d^. The horizon is 3-dimensional, with a two-dimensional space of 
generators. So long as the horizon is not itself Killing, we see that any two generators are 
related by the actions of dt and d^. 

Choose one horizon generator with affine parameter A. We can extend A to a scalar 
function on 1-L by requiring it to be invariant under dt, d^. In our case we can take A = X{w) 
since w is indeed invariant under both KVFs and is a good coordinate on Ti. 

Let k"" be the tangent to our generator associated with affine parameter A. Note that 
since dt, d^j, are also tangent to the horizon we have k ± dt, d^. We also choose any 
satisfying l"-ka = —1 and i ± dt,d(f,. We then extend k,i to vector fields defined across all of 
T-t by requiring them to be invariant under dt, d^. We then define a 'deformation tensor' Bab 
associated with fiow along the horizon generators by projecting Bab = V^fca onto the space 
orthogonal to k,i. See e.g. appendix F of 

Let us note that since dt, d^i, commute they are surface-forming, and k is orthogonal to 
this surface. So k is hypersurface orthogonal and the twist Uab = Bt^ab] vanishes. Thus 
Bab = Bba- We will use this symmetry below. 

Deviation vector fields for the horizon-generating null congruence are defined by the 
property that, when evaluated on a given horizon generator 7, they point to the same 
horizon generator 7' for all A. Let us consider a deviation vector field rj orthogonal to both 
k and i. Then (see e.g. appendix F of [53]) rj satisfies 



Since translations along dt, 9^ map one geodesic to another, both dt and d^ are deviation 
vectors. And both are orthogonal to k, i. So we may choose rji = dj for I = t,(j). Here the 
Tji are two spacetime vectors, not the components of a single vector. 

Let us now consider the set of associated inner products 



7]'Bl = kWcV"- 



(B.l) 



hij = Vi-Vj ■= VidabVj- 



(B.2) 



In any coordinate system, rj^ = dtx"" and t]'^ = d^x"". So in particular in the coordinate 
system y,w,i,(f) we have (since d^ = dt) 



hij = gif, 



(B.3) 
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i.e., this is just the induced metric on the 2-plane generated by dt, d^i, in coordinates (i, (p) or, 
equivalently for this purpose, coordinates (t, 0). So it is easy to read off from our numerics. 
But note that hu was defined to be a set of scalars, so covariant derivatives of hu are just 
coordinate derivatives. 

The evolution of hu (with respect to A, or equivalently with respect to w) is governed 
by f lB.ip . Using (]B.1|) we compute 



= 2Ba,7]'}r]'j = 2Bjj, (B.4) 

where in the last step as for fIB.Sp above we have used Bjj to denote the i, (j) components of 
Bac in the particular coordinate system y, w, i, (or equivalently the t, cj) components). 

The deformation tensor B"^^ is by construction orthogonal to k, i. Thus we may write 
I3°'^dadh = B^'^ didj. Furthermore, from fIB.Sp we have 

Bij = giaB'^'gjb = gmB'^'^gLj = hjKB'^'^hLj, (B.5) 

where K, L also range over 0, t. Thus we may safely use hjj and its inverse h^"^ to raise and 
lower indices I , J on Bjj. 

Now, the components Bjj are essentially exponentials of integrated versions of the ex- 
pansion and shear. In particular, introducing the projector Q°-f^ onto the subspace orthogonal 
to k,i (i.e., onto the space spanned by dt, d^j,) we have 

e = QabB'''' = hjjB'' = h^'Bij, (B.6) 

and 

aiJ = Bjj - -^^Ohjj, (B.7) 

where D = 4 for AdS4. 

Of course, it remains to actually find the affine parameter A used in the above definitions. 
We choose to calculate A from Raychaudhuri's equation which, in the present context, may 
be written 

JO 

— = -BabB-' - Q^'RacMk'k^. (B.8) 

As usual, the symmetries of the Riemann tensor imply that Q""^ Racbdk'^k'^ = g^'^Racbd.k'^k'^, 
which is proportional to Rcd.k^k'^. But Rab oc gab by the equations of motion, so the final 
term in (]B.8|) vanishes. Since B"''' is orthogonal to both k and i we may then write 

In terms of a general coordinate w along the generators, (1B.9P may be rearranged to yield 



^h^.hhj.j.h'j^j^h'j,j, + -^{h^'h'jj) 



X'Z{w), (B.IO) 
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where ' denotes the coordinate derivative d/dw and the last equahty defines Z{w). This 
equation is then easily solved for A in terms Z{w), which is relatively straightforward to 
extract from the numerics. 
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